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ABSTRACT 

Two recent observing campaigns provide us with moderate dispersion spectra 
of more than 230 cluster and 370 field B stars. Combining them and the spectra 
of the B stars from our previous investigations (~430 cluster and ~100 field B 
stars) yields a large, homogeneous sample for studying the rotational properties 
of B stars. We derive the projected rotational velocity Vsini, effective temper- 
ature, gravity, mass, and critical rotation speed V CT n for each star. We find that 
the average Vsini is significantly lower among field stars because they are sys- 
tematically more evolved and spun down than their cluster counterparts. The 
rotational distribution functions of V eq /V cr i t for the least evolved B stars show 
that lower mass B stars are born with a larger proportion of rapid rotators than 
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higher mass B stars. However, the upper limit of V eq /V CI i t that may separate 
normal B stars from emission line Be stars (where rotation promotes mass loss 
into a circumstellar disk) is smaller among the higher mass B stars. We compare 
the evolutionary trends of rotation (measured according to the polar gravity of 
the star) with recent models that treat internal mixing. The spin-down rates 
observed in the high mass subset (~ 9M ) agree with predictions, but the rates 
are larger for the low mass group (~ 3M ). The faster spin down in the low 
mass B stars matches well with the predictions based on conservation of angular 
momentum in individual spherical shells. Our results suggest the fastest rotators 
(that probably correspond to the emission line Be stars) are probably formed by 
evolutionary spin up (for the more massive stars) and by mass transfer in binaries 
(for the full range of B star masses). 

Subject headings: line: profiles — open clusters and associations: individual (IC 
4996, NGC 581, NGC 869, NGC 884, NGC 1893, NGC 1960, NGC 6871, NGC 
7380, NGC 7654) — stars: rotation — stars: fundamental parameters — stars: 
early-type 



1. Introduction 

The intermediate mass B-type stars are excellent targets to explore the effects of stellar 
rotation on evolution because th e consequences of rot ation are predicted to be more impor- 



tant among more massive stars (jEkstrom et al.l 120081 ) and because there are large numbers 



of B-stars to study (compar ed to the more massive O-type stars) in our region of the Galaxy. 



In several previous studies ( iHuang &: Giesll2006al lbl). we developed spectroscopic diagnostic 



tools that carefully account for the rotational distortion and gravitational darkening asso- 
ciated with rapid rotation. These studies of B stars in young clusters demonstrated that 
the stars spin down with advancing evolutionary state when the polar gravity logg po iar is 
taken as an indicator of evolutionary status. Moreover, we also noticed that the rate of 
spin down appears to be faster (relative to the main sequence lifetime) among lower mass 
stars. However, the details remain uncertain, because of the limited size of the sample and 
incomplete sample coverage. There is a pressing need for observational results on a larger 



place density ( 


Huang & Gies 
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2008 


) and internal mixing ( 
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2008; 
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2009: 
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2009 


)• 





Since then, we embarked on major spectroscopic programs in 2006 (on cluster B stars) 
and in 2008 (on field B stars). We now have much more data and much better coverage 
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of the B stars near the main sequence (MS). In this paper, we report on the new spectra 
and their analysis. We describe in §2 the two observing programs that significantly increase 
our sample size and make possible a statistical analysis of rotational trends. We outline in 
§3 how we determine the key parameters of individual stars, and we present the results in 
tabular form. In §4, we compare the projected rotational velocities of the cluster B stars 
with those of the field B stars and we argue that a systematic difference in evolutionary 
state explains the differences in their average properties. We explore the mass dependence 
of the rotational properties of young MS stars in §5, and we examine in §6 the issue of the 
minimum rotation rate necessary for rapid rotators to eject gas into a circumstellar disk and 
appear as emission line Be stars. We make a critical comparison in §7 of the observed and 
predicted changes in rotational properties with advancing evolutionary state, and we discuss 
in §8 the consequences of our work for models of the formation processes of Be stars. Our 
results are summarized in 59. 



2. Observations 

Our first run was devoted to B stars in open clusters, and we selected targets in young 
clusters in the northern sky front the WEBDA databasfi The names, ages, and number of 
B star targets of these clusters are given in Ta ble 1 . Only two of the clusters (NGC 869 and 



884) were the subject of our previous survey (IHuang fc Giesll2006al ). and the new observa- 



tions cover a different set of stars. We made the observations with the Wisconsin-Indiana- 
Yale-NOAO (WIYN) 3.5-m telescope at Kitt Peak National Observatory (KPNO) during 
2006 October 4-7 UT. We used the Hydra multi-object spectrometer to obtain moderate 
resolution spectra (R ~ 4600, using grating 1200@28.7 in second order) in the blue region 
(4000-4600A). A G7 BG-39 filter was installed to block contamination from first order flux. 
The fiber configuration for each cluster field generally consisted of 20 - 40 science stars, 5 - 
6 guide stars, and 10 - 20 random sky positions for sky background subtraction. Though we 
lost some time during the first two nights due to bad weather, we were able to obtain spectra 
of all the cluster fields on two different nights. We made a wavelength comparison exposure 
for each night and every field to secure an accurate wavelength calibration and to measure 
radial velocities with care. The errors in the radial velocity measurements (made through 
cross-correlation of the observed and synthetic line profiles) are estimated to be ~6.5 km s _1 
by comparing the two-night spectra of those cluster B stars that are not spectroscopic bi- 
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naries. After subtracting the bias level, we used the IRAFjj procedure DOHYDRA to do 
the rest of data reduction (traced and binned the spectrum, cleaned cosmic rays, calibrated 
wavelength, removed the sky background, applied the 1-D flat) to obtain the final reduced 
spectra. 

The second run focused on field B stars, i.e., those with no obvious connection with 
a known star cluster. Our sample was selected from the Smithsonian Astrophysical Obser- 
vatory (SAO) Star Catalog based on their spectral types. We avoided selecting any B star 
known to be associated with a cluster (marked as "star in cluster" in SIMBAD). We also 
checked the Deep Sky SurvejQ field images around some of our field B stars fainter than 
V — 8, and we deleted those stars found in crowded, cluster-like fields. We obtained the 
field B star spectroscopic observations with the KPNO 2.1 m telescope and Goldcam spec- 
trograph (with a 3072 x 1024 CCD detector, T3KC) during 2008 November 14 - 18 and 
2008 December 9 - 12. The spectrograph grating (G47 with 831 lines mm -1 ) was used in 
second order with a CUSO4 blocking filter, and this arrangement provided spectral coverage 
of about 900A around a central wavelength of 4400A. The slit width was set at l'/3, leading 
to a resolving power of R ~ 2400 (FWHM~1.83 A, measured in the comparison spectra). 
The integration time for each exposure was set to reach a S/N ratio > 200 in the continuum 
regions. We made comparison (HeNeAr lamp) exposures for each stellar exposure to ensure 
accurate wavelength calibration of the spectra. The errors associated with the wavelength 
calibration are approximately 6 km s -1 based upon residuals of the comparison spectra 
fit and a comparison of multi-night measurements of radial velocities of stars that are not 
spectroscopic binaries. We obtained the final reduced spectra using standard IRAF CCD 
image reduction methods (subtracted the bias level, divided the flat images, removed cosmic 
rays, and fixed the bad pixels/columns) and long-slit stellar spectrum extraction procedures 
(traced and binned the spectrum, calibrated wavelength, and normalized to the continuum). 



3. Stellar Parameters 



The procedure to d etermine the physical p arame ters of B stars is similar to that used 
in our previous studies ( jHuang fc Giesl l2006al Jbl. 120081 ) . First we measure the V sin i values 
by fitting synthetic model profiles of He I A4471 and Mg II A4481, using realistic physical 
models of rotating stars (con sidering Roche geomet ry and gravity darkening). The details 
of this step are described in iHuang fc Giesl (j2006al ). For narrow line spectra, the profile 
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fitting method may yield V sin i values that are lower than the spectral resolution. In the 
low Vsini regime, the shape change in a line profile due to rotational broadening is less 
significant. Therefore, the fitting errors increase as the V sinz values drop below the spectral 
resolution. This trend is clearly shown in the V sin i measurements of our field star sample 
(see Fig. [1]). For most of the field stars with Vsini < 40 km s _1 , the fitting errors are 
comparable to or even larger than the derived V sin i values. For the stars with V sin i > 40 
km s _1 , their V^sini values are significantly higher than the fitting errors. Thus, we use 40 
km s _1 as the lower limit of the reliable Vsini values that our profile fitting method can 
measure for our 2008 field sample. The spectral resolution of the other B star spectra we 
obtained are similar to or better than that of the 2008 field sample. Thus, the lower limit of 
reliable Vsini values for them is similar to or lower than 40 km s _1 . 

We then derive both the effective temperature (T eff ) and gravity (log g) averaged over 



the vis ible hemisphere of a rotating star by fitting the H7 profile (see details in lHuang fc Gies 



2006b|). We note that the synthesized H7 profiles were created by the spectrum synthesis 



code, SYNSPEC430, using Kurucz LTE model atmospheres. This ensures that the stellar 
parameters of the new B targets are derived in the same way as done in our previous sur- 
veys. More sophisticated non-LTE, line blanketed models yield model H7 profiles that are 
very similar t o those derived from Ku rucz LTE model atmospheres (especially for MS stars). 



For example, lLanz fc Hubenyl ( 120071 ) show that the non-LTE H7 profiles tend to be slightly 
stronger and broader than those derived from LTE atmospheres (see their Fig. 7), which 
means that our derived temperatures will slightly lower and gravities slightly higher than 
those from a non-LTE treatment. These differences are insignificant for the analyses that 
follow in subsequent sectio ns. The errors given i n our final results (see Tables 3 and 4) are 



numerical fitting errors (see lHuang fc Giesll2006al Jbl for the derivation of these errors). All the 



synthetic profiles were convolved with the instrumental broadening profiles measured from 
the corresponding comparison spectra before fitting to the observed spectra. By shifting 
the best fit profiles in wavelength to match the observed profiles, we then obtained a radial 
velocity for each night of observation. The radial velocities were transformed into the helio- 
centric frame by removing the orbital motion of the Earth (derived using the RVCORRECT 
function in IRAF). 

Stellar rotation in a close binary system is significantly affected by strong gravitational 
interactions between the components. Consequently, their rotational evolution is dominated 
by tidal effects, and we need to separate them from single stars and long period binaries 
when we study stellar rotation. We rely on the radial velocity measurements to identify close 
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binaries. If the difference in the radial velocity measurements of a star with multiple-night 
observations is larger than 13 km s _1 (~ 2a), we classify it as a spectroscopic binary (SB). We 
refer to the rest as constant velocity (probable single) stars. To check the detection efficiency 
of SBs using our 13 km s" 1 criterion, we searched The Ninth Catalogue of Spectroscopic 
Binary Orbitfl (SB9; IPourbaix et~aD l2004h for known SBs among our field B star sample. 
We found a total of 20 systems in common that are listed in Table 2, and we detected eight 
of these as SBs using the 13 km s" 1 range limit. These are mainly short orbital period SBs, 
and the detection efficiency is approximately 70% for periods P < 14 d. Our SB-detection 
efficiency quickly drops for SBs with longer orbital periods where the velocity changes over 
a 1 - 2 d interval are small. Note that although the constant velocity stars in our sample are 
sometimes called single stars in later sections, some may be SBs with longer orbital periods 
whose rotationally-related properties are probably similar to those of single stars. 

The derived gravity (logg) for a rotating star represents an average of gravity over its 
visible hemisphere. It may not be a good indicator of evolutionary status of the star because 
the effective gravity of the equatorial region is lowered by the centrifugal force induced by 
stellar rotation. Instead, we use the gravity at the poles of the star (log <? po iar) to estimate the 
age of the star since the logg po i ar value is not significantly influenced by rotation. The higher 
the Vsmi value is, the greater the difference between log g and logg po i ar . This difference, 
A log g, also depends on stellar mass, radius, and the inclination angle (i) of the spin axis. 
Ideally, if we know the inclination angle i, then we can determine A log g from the measured 
Vsmi, T e s, and logg. However, in practice, i is generally unknown. Our strategy is to 
adopt a statistical average value of A log g over the range of possible inclination angle. The 
detailed procedure is describ ed in iHuang fc Gied (j2006bl ). Based on our model simulation 
results (IHuang fc Giesll2006bl ). estimations of log g po iar are quite accurate for most situations 
(the statistical errors are 0.03 dex or less). For the extreme case where a star is spinning very 
close to the breakup velocity (K q /V^ rit > 0.95) and the inclination angle is low (i < 20°), 
the estimated logg po i ar could be off by —0.1 dex or worse, but such cases are rare. 

We can approximately deri ve the stellar mass b y assuming that the star follows a non- 
rotating star evolutionary track ( ISchaller et al.lll992[ ) in the (T e g,log^) plane with a current 
position at (T eff , log <7 po iar)- Here we assume that our derived hemisphere average T e g- is a 
reasonable choice for the whole surface average T e g. Any differences between model tracks for 
rotating and non-rotating stars will clearly influence the resultin g mass estima te. We plot the 



evolu tionary tracks for both non-rotating (ISchaller et al.lll992l ) and rotating ( lEkstrom et al. 



20081 ) models in Figure 1. This comparison shows that the differences are minimal between 
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the slow (Q/Q CT it = 0.1) and moderately fast (O/flcrit = 0.5) rotating models, while the 
evolutionary tracks of the very rapidly rotating models (Q/Q CTit = 0.9) are shifted to lower 
temperature regions occupied by the lower mass, non-rotating models. Our method may 
lead to a systematic underestimation of stellar mass by about 10% for such rapid rotators. 
Note that there is a significant shift in effective temperature between the non-rotating and 
slow rotating (f2/f2 crit = 0.1) tracks for a given mass that indicates that there are systematic 
model differences unrelated to rotation. One other error source for the estimated mass is 
related to the placement in the "zig-zag" portion of the evolutionary tracks that corresponds 
to the terminal age main sequence (TAMS; gray area in Fig. [2]). Since stars spend more time 
on the cooler portion of the track compared to the hotter, advanced stage of evolution, we 
decided to use the low temperature segment of the tracks in the shadowed area to determine 
stellar mass. The derived mass values of those more evolved stars in this region could be 
overestimated by up to 10%, although the percentage of such stars is expected to be low. 

Once we estimate stellar mass, we can derive the polar radius (r p ), and then calculate 
the critical rotational velocity using the formula for the simple Roche model (point-source 
gravitational plus rotational potential fields), 



Here we can see how the mass errors propagate into V^t errors. If mass were overestimated by 
10%, then r p would be overestimated by 5% (from logg po i ar ), and V crit would be overestimated 
by only 2.5%. Thus, the errors in mass described in the previous paragraph will result in 
relative errors in V crit of < 2.5%. 

The final derived parameters of our newly measured sample B stars are given in Table 3 
(cluster stars) and Table 4 (field stars) whose full contents are available only in the electric 
form. Table 3 lists the cluster name, stellar identification number from WEBDA, radial 
velocities, Vsini, and the derived physical parameters. The heliocentric Julian dates (HJD) 
of observation are given for each cluster in the electronic table header. Table 4 lists the star 
identifications in the SAO and Henry Draper (HD) catalogs, HJD of observation, and then 
the other quantities as in Table 3. The total numbers of B stars in our complete sample are 
summarized in Table 5 according to cluster or field membership and radial velocity variability 
status. Combining both new (first two rows in Table 5) and prior data (last two rows in 
Table 5) gives us a large and homogeneous sample (~ 1000 stars) for the statistical analyses 
presented in following sections. 




(1) 
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4. Field and Cluster Star Rotation Properties 



We argued above that the rotational properties of binary stars may not be representative 
of B-stars as a whole, and so we set aside the subsample of radial velocity variable stars 
(probably binaries). The next issue to resolve is whether B stars in the field and in clusters 
have such different rotational pro perties that the re sults for stars in the two cat egories need 
to be analyzed separately. Both lAbt et al.l ( 120021 ) and iHuang fc Gied (j2006al ) found that 
nearby field B stars appear to rotate slower on average than B stars in young open clusters. 
The authors argue that this difference results from evolutionary spin-down with time and 
the relatively greater age of the field stars. The same system a tic dif feren ce in rotatio n rates 
was a lso demonstrated in a series of papers by IStrom et al.l (120051 ) and IWolff et al.l ( 120071 . 
20081 ). These authors argue instead that the trend results from the variation in density of the 
birthplace enviro nment. To settle this issue, we recently presented an analysis of spectra of 
108 fie ld B stars (IHuang &: Giesll2008l ) that was done in the same way as our earlier cluster 
study (IHuang fc Giesll2006bl ). This work showed that both the field and cluster stars exhibit 
a similar spin-down with advanced evolution (when the stars are grouped into common 
log^p i ar bins) and that field stars are generally more evolved than the those in clusters. 
Here we repeat this comparison strategy with much larger samples. 

The fact that, on average, cluster B stars rotate faster than field B stars is clearly con- 
firmed again with our new B star samples. We show in Figure [3] the cumulative distribution 
functions of V sini for the total samples of 419 field and 524 cluster stars (all constant ve- 
locity stars). The mean V sini is 110.7 ± 4.6 km s _1 for the field sample and 145.2 ± 4.1 
km s -1 for the cluster sample. A Kolmogorov-Smirnov (K-S) test shows that there is a very 
low probability (2 x 10~ 10 ) that the two samples are drawn from the same population. 

If we plot the logg po i ar histogram for both the cluster and field B stars (Fig. |4J, we 
can clearly see the difference in evolutionary state of the two samples: the field B sample 
contains more stars with lower logg po i ar . The mean logg po i ar is 3.80 for the field B star 
sample and is 4.11 for the cluster sample. Therefore, we suggest that the slower rotation of 
our field sample is due solely to the greater influence of evolutionary spin down in a relatively 
older population. This can be demonstrated by constructing diagrams of mean Vsini as 
a function of logg po i ar for different mass groups among both the field and cluster B stars. 
We show such plots in Figure for the field (left) and cluster stars (right) for three mass 
groups: a low mass group (2 < M/M < 4; top), a moderate mass group (4 < M/M Q < 8; 
middle), and a high mass group (M/M Q > 8; bottom). For each mass group, we calculate 
the mean of Vsini for stars that are binned according to logg po i ar (over a range of 0.2 dex 
for each bin). Note that since there are only six young (logg po i ar > 4.0) field stars in the high 
mass subgroup, we consolidate two high logg po i ar bins into one as shown in Figure El If the 
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bin contains data for six or more stars, then we also calculate the standard deviation of the 
mean (shown as shaded areas in Figure [5]). We find that the evolutionary spin down trend is 
present in all subgroups of both the field and cluster samples (although with shallower slope 
among the more mas sive stars). The results are also in qualitative agreement with those of 
Abt et al.l (120021 ) and lAbtl ( 120031 ) who examined a similar number of B stars and found that 
between luminosity class V (unevolved) and class III (evolved) the mean projected rotational 
velocity declines by ~ 15% for early-B types and by w 40% for late-B types. If we compare 
the individual bins of the field stars with the corresponding ones of the cluster stars, we 
find that all of the matched pairs have similar mean V sin i (within one standard deviation) 
except for one set with 3.6 < logg po i ar < 3.8 in the mid-mass subgroup. Note, however, that 
the number of field B stars in the high mass group (> 8M ) is small (a total of only 46 
field B stars), and a larger sample would help establish the deta ils of the spin-down for this 
group. These results are very similar to our previous findings (IHuang &: Giesll2008l ) which 
were based on a smaller sample of field stars. Thus, since the spin-down trends appear to be 
similar in the field and cluster stars and since the field stars tend to be more evolved (lower 
logg po i ar ), we again conclude that the lower rotation rates of the field stars are primarily the 
result of evolutionary spin-down changes. 



In contrast to our results, IWolff et al.l (120071 ) did not detect a significant evolutionary 
change in stellar rotation among their sample of stars, and they concluded that differences 
in initial conditions and mass densi ties of star f orming regions are key to explaining the 
differences in rotational properties. IWolff et al.l (120071 ) compared rotational properties of 



B0 - B3 stars (correponding to an approximate mass range of 6 — 12M ) in young and 
old clusters and associations in both low density (p < 1M pc~ 3 ) and high density (p > 
1M Q pc~ 3 ) environments. They found that a small spin-down was evident when comparing 
the cumulative distribution functions of V^sini for young and old dense clusters (see their 
Fig. 3). This relatively small spin-down rate probably results in part because their sample 
consists mainly of more massive B-stars and because the spin-down rates may be slower (as 
a function of log g po iar) among massive B stars (see Fig. 5). Furthermore, some of the spin- 
down trend may be lost when samples are binned according to cluster age. The time scale of 
interest for evolutionary changes is the main sequence lifetim e, which depends sensitively on 
stellar mass, tms oc M~ 2 . According to lSchaller et al.l ( 119921 ). a 6M star has tms = 63 Myr 
while a 12M star has tms = 16 Myr. Thus, while the young age group of cluster stars (1 
- 6 Myr) considered by Wolff et al. consists of mainly unevolved stars, the group of older 
cluster stars (11 - 16 Myr) probably includes some evolved massive stars (69 — 100% tms for 
M = YIMq) and many unevolved, lower mass stars (17 — 25% tms for M = 6M ). Thus, 
depending on the mass distribution of stars, the rotational properties of the older group may 
show little or no evidence of evolutionary spin-down because of the mixture of evolutionary 
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states represented within the sample. 



Wolff et al.l (120071 ) found much larger differences in the rotational properties between 



stars in high and low density environments. By their criterion, all the cluster stars that we 
observed will belong to the high density group (and all the field stars to the low density 
group). We estimated cluster densities (for both the clusters listed in Table 1 and those 
observed earlier that are listed in Table 4 of Huang & Gies 2006a) from the total masses 
derived by Piskunov et al. (2008; their log Ma) and from the cluster angular radii and 
distances from Kharchenko et al. (2005; using their R cove ). We have five clusters in common 
with the sample from Wolff et al.: NGC 869, 884, 2244, 7380, and IC 1805. The density 
comparison is difficult for NGC 869 and 884 because these cl usters overlap in the sky, and 
their log Ma values are overestimated (jPiskunov et al.ll2008l ). For NGC 2244, 7380, and 
IC 1805, the derived mass densities are 5, 18, and 15 M Q pc -3 , respectively , which are 
compa rable to the corresponding results of 11, 5, and 10 M pc -3 derived by IWolff et al. 



(120071 ) using somewhat different assumptions. The mass density range of our sample of 



clusters runs from 1 to 78 M pc 3 (and higher for NGC 869 and 884), and, t h us, a ll 



the clusters we observed belong to the high density category defined by I Wolff et al.l ( 120071 ) . 
However, the distributions of V sin i in our high (cluster) and low density (field) samples show 
no evidence of systematic differences when binned according to polar gravity (see Fig. 5). 
To make as best a comparison as possible with the mass range considered by Wolff et al., 
we selected a subset of stars in the mass range 6 — 12M Q and determined the rotational 
properties as a function of polar gravity. The results (summarized in Table 6) indicate that 
there are no significant differences (within the statistical errors) between the field and cluster 
mean values. Our results do not rule out the role of environmental factors (especially among 
the more massive B-stars where our sample contains fewer targets), but given that there is 
no evidence of systematic differences between the cluster and field stars within our sample, 
we will consider the combined set of field and cluster stars together in the following sections. 



5. Young Main Sequence B Stars 

When we study stellar rotation and its evolution, there is one key aspect that we want 
to glean from observations. That is how fast stars rotate when they just start off their main 
sequence journey at the ZAMS. These initial rates are important because they show us how 
fast stars rotate at one of the most important moments in their lifetime, ZAMS, representing 
both the starting point of the relatively simple (smoothly changing) MS stage and the end 
point of the extremely complicated pre-MS stage. Accurate stellar rotation data at the 
ZAMS can help us to understand the role of angular momentum in stellar evolution both 
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before and after the ZAMS (cf. IWolff et al.ll2004f ). Theoretically, the ZAMS only defines a 



moment when hydrogen fusion is ignited in the stellar core (the temporal boundary between 
"protostar" and "star"). The evolutionary change of stellar rotation after ZAMS is expected 
to be smooth and slow. Thus, by studying the rotational properties of the least evolved B 
stars in our single (V T constant) star sample, we should be able to obtain a reliable picture 
about how fast the ZAMS stars rotate. Since the ZAMS line of theoretical stellar models in 
the mass range of B stars has basically the same surface gravity (logg ~ 4.2; see Fig. [2]), we 
decided to select all single B stars with logg po i ar > 4.15 to form a subsample that we will 
call the Young Main Sequence (YMS) stars. 

We have a total of 220 B stars in our YMS sample. The stellar rotation distribution of 
V sini/V^rit for the entire YMS sample is shown in Figure O We made a simple polynomial 
fit of the hi stogram dat a of ^sini/V^ r i t (thin line), and then we used the deconvolution algo- 



rithm from iLucyl ( 119741 ) to derive the distribution of V eq /V cr i t (thick line). Figure [6] indicates 
that B stars in our sample are born with various rotation rates. The highest probability 
density occurs around V eq /V CI it = 0.49. About 6% of newborn B stars are very slow ro- 
tators (Kq/Krit < 0.1), while about 52% of the B stars are born with a value of V eq /V crit 
between 0.4 and .8. Only 1.3% of B stars are born as very rapid rotators (Kq/Krit > 0.9). 



Wolff et al.l ( I2004j ) have already shown from their data that the specific angular momentum 
of very young B stars is spread well below the critical limit. Now we can see this more clearly 
in quantitative form in the V eq /V cr it distribution curve. 

To investigate if these statistical data depend on stellar mass, we divided our YMS 
sample into three subgroups: low mass (2 < M/M G < 4; 90 stars), middle mass (4 < 
M/Mq < 8; 92 stars), and high mass (M/M Q > 8; 38 stars). We would expect to see 
distribution curves similar to that in Figure |6] if YMS B stars in the different mass subgroups 
were born with a similar spin rate distribution. However, as shown in Figure the statistics 
clearly reveal a different story: as stellar mass increases, more and more stars are born as 
slow rotators (with lower V eq /V crit values). The fraction of slow rotators (say, V eq /V crit < 0.5) 
in each subgroup are 37%, 53%, and 84% for the low, middle, and high mass subgroup, 
respectively. 

The statistical results on the stellar rotation rates of the YMS B stars present some 
very interesting findings: 1) Even if we generally accept that most massive stars are born as 
rapid rotators, there does exist a significant fraction of newborn B stars that are very slow 
rotators (e.g., V eq /V crit < 0.1). There must exist some processes that efficiently brake these 
slow rotators early on. 2) Compared to the less massive stars, the more massive stars tend to 
be born as slow rotators. Qualitatively, two factors may contribute to this. First, we know 
that stronger stellar winds occur in more massive stars, and these can effectively remove 
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angular momentum (IMeynet fc Maederll2000f). Secondly, we know that the binary frequency 
is high among the massive O-stars (IMason et al.l 120091 ). It may be that star formation 
processes may tend to deposit angula r momentum more in orbital mot ion and less in stellar 
spin among the more massive stars (jLarsonl 120071 . |2010| ). lAbtl (120091 ). for example, found 
that stars in dense environments (usually where more massive stars are found) tend to have 
more binary systems and lower mean Vsini values. 



6. Rotational Upper Limits and Be Stars 



As discussed by iPorter &: Riviniusl ( 120031 ). the "Be" phenomenon is associated with a 
wide range of B stars whose spectra have exhibited emission features in the Balmer line 
regions at some point in time. We focus here on the so-called classical Be stars, non- 
supergiant B stars whose spectra generally show broad, photospheric, absorption lines that 
are indicative of fast rotation. To form a gaseous disk around a star and generate the emission 
features, it is generally thought that the star has to rotate very fast so that the centrifugal 
force significantly cancels the gravitational force and makes it relatively easy to launch gas 
into orbit. How fast can a B star spin before it becomes a Be star? This is a controversial topic 
in the recent literature on Be stars, and the fact that there is still no widely accepted answer 
is partially due to the difficulty of obtaining reliable V eq estimates for Be stars. Because we 
can only derive V sin i from observational data, we require measu r emen ts of a large sample 
of Be stars to deal with the sin i factor in a statistical way. lYudinl ( 120011 ) carried out a large 
rotation survey of Be stars and, by assigning V „\t to each spe ctral subtype, he obtained values 
of < V eq > I < V crit > ranging from 0.5 to 0.8. iPorterl (119961 ) investigated some Be-shell stars 
which are presumably edge-on Be stars (sinz ~ 1), and his result suggests that the Be stars, 
if not different from his Be-shell star sample, should rotate at 70~80% of V ^t- This range 



was c onfirmed by the study of Be stars in the open cluster NGC 3766 by iMcSwain et al. 



(120081 ). who compared th e observational V sin i cumulative distribution function with the 
theoretical one. However, iTownsend et al.l (120041 ) point out that very rapidly rotating stars 
may experience strong gravitational darkening towards their equatorial zones that leads to 
an underestimate of the actual Vsini value. They applied a correction for this effect to 
the data available from IChauville et al.l ( 120011 ) . and their results imply that Be stars are 
spinning close to the critical rotation rate. On t he other hand, some later studies that 



took the gravitational darkening effect into accoun t dFremat et al. 



favored relatively low Kq/Krit values for Be stars. iHowarthl ( 120071 ) points out that the high 



20051 : ICranmerl 120051 ) still 



V eq /V cr it ratio (0.95 or higher) case only requires "weak" processes to lift the surface material 
into orbit, while the low ratio case needs much more energy to send gas into orbit (the velocity 
gap between surface and orbit can reach as much as 100 km s -1 in this case). He suggests 
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that, if Be stars can form in the low ratio situation, the signatures of "strong" processes in 
stellar p hotosphere (such as "large-amplitude pulsations" ) should be easily detected. More 
recently, ICranmerl (120091 ) developed a detailed theoretical pulsational model that suggests 
that a Keplerian disk could form around Be stars with V eq /V cr i t as low as 0.6. 

Although our survey was focused primarily on the rotational properties of normal B 
stars, we can use our sample to investigate this issue by considering the upper limit of 
Kq/Krit for those stars that probably are rotating at rates slightly below that required to 
form a Be disk. To achieve this goal, we first removed all the Be stars from our single 
(V x constant) star sample by excluding already known Be stars and those with Balmer 
emission or deep, narrow, absorption lines (a signature of Be-shell stars) in our spectra. 
This process led to a culling of 61 Be stars, leaving a total of 894 stars in the remaining 
sample (presumably all non-emission line B stars). These Be stars are identified in the last 
column of Tables 3 and 4. Note that our approach here makes the tacit assumption that the 
processes leading to disk formation require a certain minimum value of the ratio of equatorial 
velocity to critical rotation. Ideally, we could test this idea by finding a minimum value of 
this ratio among the Be stars, but unfortunately, it is difficult to apply our methods to Be 
stars because of emission contamination in the profiles and the small number of Be stars we 
observed. Instead, we will simply compare below our results on the fastest rotating, non- 
emission line stars with results on the Be stars from other investigations. Since the process 
of disk formation p robably requires rapid rotation and some other mechanical energy source 
(possibly pulsation; iRivinius et al.ll2003l ; ICranmerll2009l ) and since the properties of the latter 
may vary from star to star, we caution that the V eq /V cv \ t criterion for disk formation will 
probably apply only in some average sense. 

We can approach the problem of finding a nominal upper limit on V eq /V CT i t in two ways. 
First, we can consider the observed upper limit in the distribution of V smi/V crit . One could 
simply adopt the largest derived value of this ratio as the upper limit for V eq /V CTit , but because 
the errors associated with any given estimate are not negligible, we prefer to form a mean 
from a subset of stars with the largest observed ratios. Clearly, the larger the sample, the 
more well defined will be the upper limit boundary of the distribution. Those stars that help 
us set the upper limit will have a Kq/Krit ratio close to the limit and will have an inclination 
angle close to 90° (i.e., sinz ^ 1). If we choose the top 4% of stars that have the highest 
V eq /V cr it values (the fastest rotators) and assume that their spin axes are randomly pointed 
in space, then about 1/4 of these stars will have sini > 0.97. Thus, we expect that those 
stars in the top 1% of the V sini/V CTit distribution probably correspond to very fast rotators 
(in the top 4%) with an edge-on orientation (sini > 0.97). Therefore, for our sample of about 
900 stars, the mean ratio from the top nine fast rotators of V smi/V CV i t = 0.91 provides an 
initial estimate of the limit separating the non-emission line and emission-line stars. We 
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list in the lower row of Table 7 this estimate plus the unique maximum and 3% subsample 
estimates for the upper limit from all the ^sinz/V^ r i t measurements. 



A second approach is to use the deconvolution method of iLucyl (119741 ) to derive the 
Kq/Krit distribution from the V sini/V CV i t distribution (as shown earlier in Fig. |6]and[7|), 
and then determine the value of V eq /V cr \t that marks the boundary for a given percentage 
of the distribution. The final row of Table 7 gives these boundaries for the top 4% and 
0.2% of the reconstructed Kq/Krit distribution for the full sample of non-emission line stars. 
The fact that the 4% boundary ratio (Kq/Krit = 0.88) is quite similar to the estimate from 
the top 1% of the V sxni/V CT it distribution (< V sini/V^rit >i%= 0.91) indicates that both 
methods leads to similar and presumably reliable results. 

However, we found in §5 that the rotational properties of the youngest B stars in our 
sample are strongly dependent upon the mass subgroup (slower among the more massive 
stars; see Fig. [7]), and thus, it is important to consider the upper limit of V eq /V CT i t among 
samples of comparable mass. In order to obtain meaningful results from each mass subsam- 
ple, we need to keep the population large enough so that 1% of their content comprises at 
least 2 stars. Thus, we divided our total of 894 stars into four subsamples based upon mass. 
The mass range of each subsample is given in the first column of Table 7. We made the 
same two statistical estimates for the upper limit of V eq /V cr i t for each subsample. All the 
key statistical results are given in columns (2) to (7) of Table 7. The estimates based on 
the maximum value (from one star; column 3) and the mean value of < V smi/V CI i t > for 
the top 1% of the stars (based on 2-3 stars; column 4) suffer from small number statistical 
errors. Instead, estimates from < V smi/V cr it > of the top 3% of stars (5-8 stars; column 5) 
and from the boundary of V eq /V cr i t for the top 4% of the rapid rotators (6-10 stars; column 
6) provide more reliable results. 

We illustrate the rotational properties of the fastest spinning, non-emission line stars in 
two sets of panels in Figure[8l The left hand panels show the distributions in the (V sin i/V CT it, 
log g po iar) plane with a dotted line marking the mean < ^sinz/V cr i t > of the top 3% of stars, 
while the right hand panels show the histograms of V sin i/V CI i t and the deconvolved functions 
Kq/Krit (format similar to Fig. Eland [7]). The dotted and dashed lines in the right hand 
panels indicate the boundaries in the Kq/Krit distribution for the top 4% and 0.2% of the 
sample, respectively. These panels are arranged in four rows corresponding to the four mass 
subgroups (from low mass at the top to high mass at the bottom). 

We can immediately see from Figure [S] that the rotational upper limit for non-emission 
line B stars does indeed depend on stellar mass. Estimates of the upper limit from both 
< V sini/Krit >3% and (Kq/T4rit)4% show a consistent drop from about 0.92 in the 3M 
subsample to about 0.56 in the 9M & group. In principle, we might take the rotational upper 



- 15 - 



limit as the value of Kq/Krit where the deconvolved function reaches zero probability in 
Figure However, the detailed shape of this function at the top end depends critically on 
the small numbers in the last few populated bins of V sin i/Krit- Thus, we arbitrarily adopted 
upper limit as the boundary of V eq /V crit for the top 0.2% of the deconvolved distribution, 
since this value is somewhat better constrained and corresponds to the case where we might 
find or 1 star in our subsamples. It is hard to estimate the error bars of the adopted 
upper limit of V eq /V CTit from regular methods. However, since the < V smi/V crit > 3 % and 
(K q /Krit)4% statistical values in Table 7 are quite reliable and comparable, their mean can 
can serve as the lower boundary to the upper limit of K q /V^ r i t . Then a working estimate of 
the error in the upper limit can be formed by taking the difference between (Kq/Krit)o.2% 
and the mean of < V sini/Krit >3% an d (V eq /V a it)4%. Our final estimates of the rotational 
upper limits are given in column 7 of Table 7. 

These results suggest that the rotational threshold required to create a Be star disk 
varies dramatically through the B star mass range. For example, we find that low mass B 
stars (M < 4M Q , or spectral subtype later than B6 for MS stars) may rotate extremely 
fast (Ksq/Krit ~ 0.96) without creating an outflowing disk. This may explain wh y the 
fast rotator, Regulus (below the threshold at Kq/Krit = 0.86; iMcAlister et al.ll2005l ). has 
never exhibited any Balmer line emission. For stellar masses > 8.6M (spectral subtype 
B2 or earlier for MS stars), the threshold drops to V eq /V cr i t ~ 0.63. This result may help 
resolve the controversy surro unding the massive Be star, Achernar (~ 7M Q ). Interferometric 
observations of Achernar by iDomiciano de Souza et al.l (120031 ) showed that the star has an 
exceptionally large rotational oblateness with R eq /R vo \ai = 1-56, which is above the value 
ex pected (1.5) for a star r otating at the critical rate. However, a recent spectroscopic analysis 
by IVinicius et al.l (120061 ) indicates that Ache rnar may rotate at 0/f^ cr jt ~ 0.8, which is 
equivalent to V cq /V cr i t ~ 0.65 (cf. Fig. 10 in lEkstrom et al.l 120081 ) . Note that this rate is 
probably near to the implied threshold we find for disk formation among the massive B 
stars, so Achernar's rotation rate and status as a Be star are probably consistent with our 
results. 



The ma ss dependence of the rotation limit was generally ignored in past with the notable 
exception of ICranmerl (120051 ) who reported a very similar V eq /V cr n dependence on spectral- 
type. He applied Monte Carlo forward modeling to a sample of 462 classical Be stars to test 
various V eq /V cri t distribution functions against the observed V smi/V CT i t distribution. His 
conclusion was that hot B stars (T e g > 21000K) tend to have a low threshold of V eq /V CT i t 
("well below unity") while cool B stars (T c g < 21000K) have an increas ing threshold o f 
Kq/Krit (up to unity) as "T eff decreases to the end of the B spectral range". iHowarthl (120071 ) 
raised two concerns about these results: 1) they were based on very heterogeneous data 
sources, and 2) V CT it was assigned to individual stars according to spectral type. However, 
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these concerns do not apply to our sample, and the overall good agreement between our 
result (based on normal B stars) and Cranmer's (based on classical Be stars) convinces us 
that the rotational limit for Be star formation is indeed mass dependent. 

We should note that most of the B stars in our sample were retained for this analysis 
based upon an examination of only a few spectra for B aimer line emission. However, Be stars 
are intrinsically variable in nature (IPorter fc Riviniud 120031 ). and in some cases ( "tran sient 
Be stars") the disk emission may even disappear at times (cf. iMcSwain et al.l 120091 ) . It 
is possible that some of the fast rotators used in our analysis are transient Be stars that 
we observed during an emission-free epoch. If such Be stars do exist among those high 
^sini/\4 r it stars that we used to determine the rotational upper limit, then we may have 
overestimated the limits. We list in Table 8 all the stars with high ^sinz/V^ r i t that helped 
define our derived upper limits, and we recommend spectroscopic Ha monitoring of these 
targets to verify their status as normal (non-emission line) B stars. 



7. Rotational Evolution of Main Sequence B Stars 



Theoretical studies (IHeger &: Langerl l2000t iMeynet fe Maederl I200CH ) demonstrate that 



the evolutionary paths of rapidly rotating, massive stars can be very different from those 
for non-rotating stars. Rapid rotation can trigger strong interior mixing, extend the core 
hydrogen-burning lifetime, significantly alter the luminosity, and change the chemical com- 
position of the o uter layers over t ime. The most detailed models to date are presented in a 



seminal paper by lEkstrom et al.l ( 120081 ) (=EMMB), who show how the equatorial velocities 
and surface abundances change with time according to the initial rotation rate, metallicity, 
and mass. We have already shown in §4 (Fig. E]) that the mean rotational properties do 
vary with logg po i ar , our measurement of evolutionary state. In this section, we re-sample our 
rotational data into two mass categories from the EMMB model grid, and we compare the 
observed changes with log g po iar to predictions from the EMMB models that are based upon 
an initial rotational distribution derived from the youngest MS stars. 

Our goal is to compare our rotational results with the EMMB models for solar metallicity 
(a reasonable assumption for nearby, Galactic, B stars). The EMMB models are presented 
for masses of 3, 9, 20, and 60M Q , and since massive stars with M > 15M Q are sparse in our 
sample, we are limited to making the comparison with the EMMB models for 3M Q and 9M . 
We selected stars with 2 < M/M Q < 4 (262 stars; < M >= 3.2M ) for the 3M group and 
those with 7 < M/M & < 13 (241 stars; < M >= 9.2M ) for the 9M group. All the known 
and suspected Be stars were omitted from these samples because we cannot derive reliable 
T e fj and logg parameters from the emission-blended H7 profile. For each group, we first 
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derived the rotational distribution curve (probability density vs. V eq /V cr i t ) of the YMS stars 
(logg po i ar > 4.15; 90 stars in the 3 M group, and 49 stars in the 9 M group) using the 
same method described in §5. These functions describe the initial distribution of Kq/Knt- 
Then the EMMB theoretical models can be used to show how these distributions should 
change with time and, consequently, how < V smi/V cr i t > should change with time. This 
mean ratio can be directly compared to our observational results as a function of logg po i ar 
to investigate evolutionary changes. 

The EMMB models provide the key data on how stellar rotation evolves with time for 
initial (ZAMS) values of angular rotational rate f2/f2 cri t = 0.1, 0.3, 0.5, 0.7, 0.8, 0.9, and 
0.99. Figure [9] shows how the V eq /V CT i t ratio changes with logg po iar in the models for each 
of these starting rotational rates. We see that the EMMB models predict that a sharp drop 
occurs at the outset for stars with higher initial V eq /V crit values. This results from a rapid 
redistribution of angular momentum as the models transform from an assumed solid body 
rotation state at ZAMS into a differentially rotating configuration. Our YMS star sample 
is best compared to the models that have stabilized and that are characterized by a mean 
logg po i ar similar to that of our YMS samples. Thus, we adopted logg po i ar = 4.22 as the 
nominal starting point for the model calculations. Then the deconvolved distribution of 
Kq/Krit derived from the YMS stars maps into a probability density for each of the EMMB 
tracks shown in Figure [9) Note, however, that at this stage of evolution (logg po iar = 4.22), 
the EMMB tracks only sample a range from V eq /V cr it = 0.0 to 0.7, and we need to extrapolate 
beyond these tracks for the portion of the rotational distribution with the fastest spins. 

We then performed a numerical integration in small steps of A log g po iar to find how 
the probability becomes redistributed in V eq /V cv \t as evolution advances. From inspection 
of Figure [9j the part of the distribution for the slowly spinning stars will be shifted only 
slightly higher, while the porti ons of the distribution for the fastest rotators will be shifted 



to significantly higher V eq /V CTit . lEkstrom et al.l (120081 ) argue that these rapidly rotating stars 



will evolve towards the critical limit, where they may become Be stars and drain their angular 
momentum into a circumstellar disk. However, we purposely removed the Be stars from 
our observational sample, so we need to truncate the fastest portion of the evolved model 
distributions of V eq /V CT i t in order to make a valid comparison with the observed distributions. 
From our results in §6 (see Table 7), we estimate that the upper limit boundary separating 
the non-emission line and Be stars occurs near V eq /V cr it = 0.96 for a 3M Q star and near 
V eq /V cr it = 0.75 for a 9M star. These limits are indicated by dotted, horizontal lines in 
Figure |HJ Thus, we treated any part of the model Kq/Krit distribution that was transformed 
above these limits as removed from the sample. We then determined the mean value of 
< V eq /V CT it > from the remaining model distribution at each increment of logg po i ar (over 
the range 4.22 > logg po iar > 3.5; indicated by vertical dotted lines in Fig. [9]). Finally, this 
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statistic was multiplied by 7r/4 (assuming random orientation) to obtain a model prediction 
of < V smi/V CTit >. These calculated, mean rotational rates are plotted as thick solid lines 
in Figure [TU] for the 3M and 9M Q cases, and the predictions suggest that the sample means 
should be almost constant over the main sequence lifetime. This may appear somewhat 
surprising, given the general trend of increasing < V eq /V cr i t > for the tracks shown in Figure|9j 
but recall that the rapid rotators are systematically removed from the sample distribution 
once they exceed the nominal Be star boundary. This removal of the high end part of the 
distribution tends to cancel the effect of the small increase in < Kq/Krit > f° r the lower 
velocity portion and results in almost no change in the distribution mean value. 

We also show in Figure ITUl the observed distributions of ^sini/V^ r it for the two mass 
samples. The mean values of V^sini/V^ rit of the observed stars in each 0.2 dex bin of log g po i a r 
are plotted as horizontal bars, and their associated standard deviations of the mean are 
shown as gray shaded regions (as in Fig. [5]). The differences between the two mass groups 
are striking. The lower mass B stars (top panel of Fig. [TO]) rotate much faster than do 
the higher mass B stars (lower panel of Fig. ITU]) at the beginning of the MS stage, but 
the lower mass stars spin down much "faster" with respect to logg po i ar than do the higher 
mass stars. At the end of the MS stage, both low and high mass B stars have r oughly the 



same < V smi/V rrit > ■ Our results are in qualitative agreement with those of lAbt et al. 



( ]2002l ) and lAbtl (j2003f ) who examined a similar number of B stars and found that between 
luminosity class V (unevolved) and class III (evolved) the mean projected rotational velocity 
< V sin i > changes little for early-B types and drops significantly for late-B types. 

A visual appraisal of the predicted and observed evolutionary trends indicates that the 
EMMB model prediction makes a good match of the observational results for the high mass, 
B star group but fails to reproduce the observed spin-down of the lower mass group. The 
evolutionary trends in the EMMB models are directly related to their formulation of the an- 
gular momentum redistribution in the stellar interior caused by m ixing. The EMMB model s 



occupy the middle ground between two extreme cases discussed by lOke fc Greensteinl ( 11954J ) . 



In case A of lOke &: Greensteinl ( 119541 ) . the angular momentum is fully fixed mixed within 
the interior and the star rotates as a solid body. As the star evolves with time, it becomes 
more concentrated towards its center (and its moment of inertia gradually decreases), and 
consequently the outer layers spin up in terms of V eq /V CTit . Case B, on the other hand, imag- 
ines the star as a sequence of spherical mass shells in which each shell maintains its angular 
momentum. Then the surface rotational velocity will drop in proportion to the increase in 
radius with advanced evolution. Case A models, with their greater spin up, would tend 
to increase faster with advancing evolution than the EMMB models shown in Figure [TUJ 
while case B models, with more spin down, would decrease faster with advancing age. Since 
the observations suggest more spin down for the 3M group than predicted by the EMMB 
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model, it is worthwhile considering the case B alternative. 

In case B, we suppose that angular momentum is conserved in the individual spherical 
shells of a rotating star. Then, conservation of angular momentum in the outermost layer 
(where we measure V sini) directly leads to V cq oc r _1 . Since V cr i t oc r _1//2 , the ratio will vary 
as V eq /V CV it oc r -1 / 2 oc g l J^ &r The predicted trends of spin-down with logg po iar for the case B 
model are plotted as dashed lines in Figure [9j Then we can calculate how the distributions of 
Kq/Krit will change with advanced evolution according to the case B model and determine 
the mean < V svni/V cr - lt > in the same way as we did for the EMMB models. The case B 
model predictions are plotted as thick, dotted lines in Figure [TOj 

Surprisingly, the prediction from the simple case B model makes a good fit of the spin- 
down trend observed for the low mass group. For the high mass group, the fits from the 
case B and EMMB models are both reasonably good (except for the lowest gravity bin at 
log <7 P oiar ~ 3.5 where the EMMB model fit is better). The small difference in the predictions 
between the EMMB and case B models is due to the placement of most of the high mass 
stars in the lower V eq /V CTit region (bottom inset panel of Fig. [TD]) where the evolutionary 
changes of both models are modest (see Fig. Note that the case B model ignores any 
angular momentum loss from the surface caused by stellar winds. Such loss is negligible 
for the weak winds of the lower mass stars, but inclusion of the influence of the stronger 
winds of the higher mass stars would cause a somewhat larger spin-down than illustrated in 
Figure [TU] and might result in a more discrepant fit. 

We suspect that the failure of the EMMB models to match the faster spin-down we 
observe in the lower mass stars is related to missing physical processes that involve mag- 
netism. The transport of internal angular moment um is a complex process that iny olves 



mixing, magnetic fields , and p ossibly gravity waves (jTalonll2008l ; iDecressin et al.ll2009l ). For 



example, iMaeder et al.l ( 120091 ) describe how differential rotation, meridional circulation, and 
mixing can generate a magnetic dynamo in some circumstances. The magnetic braking as- 
sociated with the magnetic field and mass loss may then lead to a significant spin down of 
the star. We speculate that such magnetic braking is more apparent in the lower mass stars 
because their MS lifetimes are relatively much longer and consequently magnetic braking 
has a greater cumulative effect. 



8. Consequences for Be Star Formation 



Our results offer some insight about the evolutionary stages where we might expect 
to find the rapid rotation associated with disk formation in Be stars. There are three 
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broad explanations for the origin of the rapid rotating Be stars: 1) they were born as rapid 
rotators; 2) angular momentu m redistribution with evolution causes their outer layers to 
spin up ( lEkstrom et al.l 120081 ); and 3) they are the mass gainers th at were spun up by a 
past episode of Roche lobe overflow in an interacting binary system (IPols et al.lll99ll ). We 
showed in §5 that most young B stars have rotational velocities that are well below the limit 
for Be star formation (compare Fig. [7] for the YMS stars with Fig. [8] for the upper limits 
of non-emission line stars). Consequently, we think that only a small fraction of Be stars 
were born as rapid rotators, and below we will confront the expectations from the other two 
explanations with our results. 

We showed in the last section (see Fig. [TOj lower panel) that the EMMB model pre- 
dictions about evolutionary changes in < V sinz/V^ rit > are consistent with the observed 
means for the higher mass stars. The EMMB models suggest that many of the stars in the 
upper part of the rotational distribution will evolve towards crit ical rotation (Fig. [9j lower 
panel) and presumably become Be stars. iMcSwain fc Gied (120051 ) found that the Be fraction 
increases as high mass B stars evolve. Their data show that t he Be fraction of e arly B stars 



is 4.7% while that of the evolved B stars increases to 7.4%. IZorec et al.l (120051 ) also found 
that the Be star fraction increases be tween the early to late MS stages in Galactic stars in 
the mass range from 5 to 12M , and iMartayan et al.l ( 120061 ) confirmed a similar trend for 
Be stars in the Large Magellanic Cloud (LMC). All these results are consistent with the 
EMMB picture that rapid rotation develops through the angular momentum redistribution 
that occurs as the stars approach the TAMS. 

Using the EMMB model shown in the bottom panel of Figure [9] and our YMS data 
in the embedded portion of Figure [TDl we can calculate how the Be fraction might change 
with time. According to the EMMB models, Kq/Krit will generally increase with decreasing 
logg po i ar , and once a star reaches the critical limit, it will remain there until the completion 
of core H burning. Suppose that a fraction x of the sample are Be stars at the outset of the 
YMS star distribution (i.e., born as rapid rotators). Then, we can use the EMMB tracks 
to determine the number of high mass stars that should cross the Be formation limit as a 
function of evolutionary logg po i ar . The calculation indicates that at logg po i ar = 4.1, the Be 
fraction is only about x + 3% (formed from all the YMS B stars with V eq /V CT i t > 0.7, see 
the embedded plot of Fig. (TD]); at logg po i ar = 3.9, the Be fraction increases to x + 12% (from 
the YMS B stars with V eq /V cr it > 0.6); and at logg po i ar = 3.6, the Be fraction increases to 
x + 22% (from the YMS B star s with Kq/Krit > 0.5). This final estimate is in reasonable 
agreement with the results from IMcSwain fc Gied (120051 ) (Be fraction for ev olved B stars of 
7.4% ), given that the Be detection efficiency may be as low as 25 - 50% (IMcSwain et al. 
2008h . 
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However, the situation is different for the lower mass B stars. We found that the mean 
spin rates of the lower mass stars decrease faster than predicted by the EMMB model (Fig. [TO} 
upper panel), which probably indicates that these stars evolve to lower Kq/Krit (dashed lines 
in Fig. EH upper panel). This conclusion is supported by the statistics of late type Be stars. 
If the EMMB models did apply to the lower mass B stars, then we can calculate the Be star 
fraction with time in the same way as we did above for the higher mass stars. Since the 
lower mass star distribution has a higher mean velocity, a larger proportion of the sample 
are predicted to end up as Be stars (see Fig. [HJ upper panel). For example, all the low mass 
stars with V eq /V CT i t > 0.6 at \ogg po \ ar = 4.22 will be transformed into Be stars during the 
MS stage even though the boundary for Be formation is very high. According to the top 
panel of Fig. [7J about 45% of low mass YMS stars have V eq /V CT i t > 0.6, a nd all of them are 



predic ted to become Be stars before they reach the end of the MS stage. iMcSwain &: Gies 



( 120051 ) found that the Be fraction of the late B stars is only about 0.9%, which is much lower 
than the prediction from the EMMB model. Thus, the low observed Be fraction among the 
lower mass stars offers indirect evidence that these stars spin down faster than predicted. 

This line of reasoning suggests that lower mass stars spin down with time. If this is also 
true of low mass Be stars, then any such star must be in a state where it obtained its high 
spin very recently, i.e., it is one of the rare newborn stars with fast spin or it was spun up 
through binary mass transfer. We note, for example, t hat Regulus was r ecently discovered to 



be a spectroscopic binary with a low mass co mpanion (IGies et al.ll2008l) that is probably the 



white dwarf remnant of a prior mass donor (IRappaport et all 120091 ). We suspect that this 



past mass transfer did spin up Regulus to the critical rate and that Regulus may have been 
a Be star in the past. However, with the passage of time, it has now spun down to a level 
below the critical level and its rotational rate decline will likely continue. If this scenario 
for Regulus is correct, then we suspect that most of the low mass Be stars and other rapid 
rotators were formed through mass exchange. Their binary character may be tested through 
dedicated spectroscopy and radial velocity measurements. 



9. Conclusions 

Two recent observing campaigns on cluster and field B star provide us with new rota- 
tional data on about 600 targets. Combined with data obtained in our previous surveys, we 
now have a very large sample of B stars that is ideal to investigate many questions about 
the rotation of massive stars. Following the analysis methods we used in previous study, we 
derived Vsmi, T e g, and logg po i ar (considering rotational distortion and gravity darkening) 
and estimated the stellar mass and V CI it for each star. We summarize our findings below: 
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1) The radial velocity measurements allow us to identify the short period (P < 14 d) 
spectroscopic binaries (SB). The fraction of SBs in our field and cluster samples are similar. 
Considering the efficiency of detection, we estimate that about 19% of B stars are short 
perio d SBs. This result is consiste nt with many previous studies of binary frequency of B 



stars (IWolff Hl978l : lAbt et al.lll990h . 



2) The mean Vsini for the field sample (441 stars) is slower than that for the cluster 
sample (557 stars), confirming results from previous studies. By comparing stars with similar 
evolutionary status, we find that the stars in these two samples have similar rotational 
properties when plotted as a function of logg po i a r- The mean rotation rate of the field 
B sample is slower than that for the cluster sample because the former contains a larger 
fraction of evolved stars. 

3) The rotation distribution curves based on young stars with logg po i ar > 4.15 suggest 
that massive stars are born at various rotation rates, including some very slow rotators 
(Kq/Krit < 0.1). The mass dependence suggests that higher mass B stars may preferentially 
experience angular momentum loss processes during and after formation. The low mass 
stars are born with more rapid rotators than high mass stars if the stellar rotation rate is 
evaluated by V eq /V CT it. 

4) The statistics based on the normal B (non-Be) stars in our sample indicates that 
low mass B stars (< 4M ) may require a high threshold of Kq/Krit > 0.96 to become Be 
stars. As stellar mass increases, this threshold decreases, dropping to 0.64 for B stars with 
M > 8.6M . This implies that the mass loss processes leading to disk formation may be 
very different for low and high mass Be stars. 

5) Comparing; w ith modern evolutionary models of rotating stars (for 3 and 9M & from 



Ekstrom et al.l 120081 ) . the apparent evolutionary trends of < V smi/V crit > are in good 
agreement for the high mass B stars, but the data for low mass B stars shows a more 
pronounced spin-down trend than predicted. 

6) Predictions for the fractions of rapid rotators and Be stars produced by the redistri- 
bution of angular momentum with evolution agree with observations for the higher mass B 
stars but vastly overestimate the Be population for the lower mass stars. The greater than 
expected spin-down of the lower mass stars explains this discrepancy and suggests that most 
of the low mass Be stars were probably spun up recently. 
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Table 1. Open Clusters Observed 





log Age 


Number of 


Name 


(y) 


Stars Observed 


IC 4996 


6.95 


16 


NGC 581 


7.34 


22 


NGC 869 


7.07 


35 


NGC 884 


7.03 


34 


NGC 1893 


7.03 


22 


NGC 1960 


7.47 


27 


NGC 6871 


6.96 


25 


NGC 7380 


7.08 


22 


NGC 7654 


7.76 


31 



Table 2. Spectroscopic Binaries Common to the Field B Star Sample and SB9 



SB9 Our Measurements 

SB9 P K{ ~K 2 7 K-(N1) V r (N2) |AV^| 

Index HD (d) e (km s _1 ) (km s _1 ) (km s _1 ) (km s^ 1 ) (km s _1 ) (km s _1 ) Detected 



131 


16219 


2.1 


0.05 


21 




12 


4 


28 


24 


y 


1353 


209961 


2.2 


0.03 


122 




— 13 


— 14 


— 19 


5 


n 


189 


23466 


2.4 


0.0 


23 




1 ( 


—7 


33 


39 


y 


484 


65041 


2.8 


0.30 


34 




— 13 


26 


—29 


55 


y 


1423 


O 1 O A AT 

218407 


3.3 


U.z4 


86 




1 r\ 

— 10 


— 2 


—69 


e +7 
D I 


y 


146 


17543 


3.9 


0.04 


25 




o 
b 


7 


— 11 


1 1 


y 


313 


34762 


5.4 


0.08 


27 




6 


-1 


-17 


16 


V 


1424 


218440 


7.3 


0.38 


88 


147 


-5 


-22 


-10 


13 


11 


394 


44172 


8.2 


0.0 


25 


48 


-6 


-24 


-17 


7 


n 


1409 


216916 


12.1 


0.05 


24 




-12 


-12 


3 


15 


y 


375 


41040 


14.6 


0.39 


35 


39 


13 


11 


9 


2 


n 


17 


1976 


25.4 


0.12 


24 




-10 


12 


2 


10 


n 


44 


4382 


33.8 


0.41 


16 




-4 


5 


-8 


13 


y 


2401 


20315 


36.5 


0.3 


20 




4 


-9 


-9 





n 


297 


32990 


58.3 


0.19 


37 




16 


36 


34 


2 


n 


19 


2054 


48.3 


0.38 


30 




3 


-14 


-7 


7 


n 


95 


11529 


69.9 


0.30 


30 




-25 


7 


4 


3 


n 


29 


3322 


399.6 


0.57 


8 




4 


11 


8 


3 


n 


66 


7374 


800.9 


0.31 


4 




-15 


-21 


-16 


4 


n 


36 


3901 


940.2 


0.4 


12 




-11 


-4 


-2 


2 


n 



The Ninth Catalogue of Spectroscopic Binary Orbits ([Pourbaix et al. 



2004). 



Table 3. Stellar Parameters of New Cluster B Stars 



Cluster WEBDA V r (Nl) a V r (N2) a Vsini AVsini T cff AT eff Mass V crit 

Name ID (km s _1 ) (km s _1 ) (km s _1 ) (km s _1 ) (K) (K) logg Alogg log3 po i ar (M e ) (km s _1 ) Comment 

IC4996 4 -6 -6 25 12 10508 50 4.183 0.022 4.193 2.6 392 

IC4996 6 -20 -26 156 7 10270 50 3.903 0.010 4.032 2.7 361 

IC4996 8 -5 -13 64 9 26014 550 4.262 0.062 4.273 10.4 582 



a The HJD times for V r are given in the online version of this tabic. 



Table 4. Stellar Parameters of New Field B Stars 



HJD(Nl) Vr(Nl) HJD(N2) V r (N2) Vsini AVsini T cB AT cff Mass V clit 

SAO HD 2454000+ (km s" 1 ) 2454000+ (km s" 1 ) (km s" 1 ) (km s" 1 ) (K) (K) logg Alogg logc/ polar (M ) (km s" 1 ) Comment 



4080 1359 790.669 1.5 791.659 -2.3 203 12 12000 50 3.834 0.012 4.007 3.4 376 

4165 3366 790.671 2.8 791.662 -5.1 3 21 17729 250 3.886 0.027 3.887 6.3 411 

4226 4382 790.675 5.0 791.664 -8.3 72 31 12612 150 3.299 0.029 3.393 4.9 290 



I 

CO 

o 

I 



31 



Table 5. Summary of the B Star Samples 



Sample 


N(V T const.) 


N(V T var.) 


iV(Total) 


New Field 


332 


43 


375 


New Cluster 


205 


27 


234 a 


Prior Field b 


91 


13 


104 


Prior Cluster 


327 


78 


433 a 



a Includes some cases where the V T variation status is 
unknown since only one measurement is available. 



iHuang & Giesl (120081) 



Huang k Giesl (12006bl ) 
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Table 6. Rotational Data for the Field and Cluster B Stars with 6 < M/M & < 12 







Field 






Cluster 




Range of 


No. of 


< V sin i > 




No. of 


< V sa\i > 




log ffpolar 


Stars 


(km s _1 ) 


(km s _1 ) 


Stars 


(km s _1 ) 


(km s _1 ) 


3.2-3.4 


5 


26.2 


25.8 


1 


0.0 




3.4-3.6 


17 


76.5 


18.6 


10 


102.1 


27.9 


3.6-3.8 


19 


77.4 


19.8 


17 


95.5 


17.6 


3.8-4.0 


18 


107.7 


26.5 


82 


115.1 


9.6 


4.0-4.4 


15 


163.3 


34.5 


136 


135.1 


9.1 



Table 7. Statistical Results for the Fastest Rotating Normal B Stars 



Mass Range 




From V smi/V CI i t distribution 


From Kq/Vcrit 


distribution (V cq /V cri t) uppcr limit 




N 


(Vsinz/V„it)ma X < V Sm ij V CI it > 1% < V Sini/Vcrit > 3 % 


(V r oq /V r crit)4% 


(V C q/V cr it)o 2% for non-Be stars 


(1) 


(2) 


(3) (4) (5) 


(6) 


(7) (8) 



2.2-4.0 255 

4.0-5.4 225 

5.4-8.6 255 

8.6— max 159 

all 894 



0.98 
0.91 
0.82 
0.59 
0.98 



0.96 
0.88 
0.81 
0.58 
0.91 



0.91 
0.84 
0.76 
0.55 
0.85 



0.93 
0.87 
0.79 
0.57 
0.88 



0.99 
0.94 
0.86 
0.63 
0.95 



0.99 
0.94 



+0.01 
-0.07 
+0.06 



-0.08 

0.86j-° 08 



0.63 



-0.08 
+0.07 
-0.07 



I 

CO 
CO 

I 
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Id Die O. 


The Highest V smi/V crit 


Stctrs 


Star a 


Mass (M Q ) 


V sin i lVrv\t 


log ffpolar 


HD236849 


3.6 


0.98 


3.652 


NGC7160-940 


3.8 


0.96 


3.714 


HD26793 


3.4 


0.93 


3.965 


HD12352 


2.9 


0.92 


3.976 


HD8053 


3.5 


0.91 


3.942 


HD9766 


3.9 


0.89 


3.364 


HD87901 


3.5 


0.87 


3.950 


HD53744 


3.3 


0.86 


3.818 


HD52860 


4.0 


0.91 


3.380 


NGC869-1417 


4.3 


0.86 


4.175 


NGC7654-1174 


4.2 


0.85 


3.904 


NGC7654-1272 


4.1 


0.84 


3.949 


HD5882 


4.6 


0.82 


4.131 


HD2626 


4.3 


0.80 


3.753 


NGC1960-91 


4.2 


0.79 


4.195 


NGC1893-256 


5.5 


0.82 


4.267 


NGC6871-14 


6.2 


0.81 


4.170 


NGC3293-28 


8.3 


0.79 


3.933 


NGC4755-139 


5.6 


0.75 


4.266 


NGC3293-29 


7.2 


0.74 


4.108 


NGC6871-16 


6.1 


0.74 


4.151 


NGC884-2255 


7.8 


0.73 


3.901 


NGC884-2622 


6.3 


0.73 


4.212 


NGC884-2190 


9.8 


0.59 


4.102 


HD52918 


10.1 


0.57 


3.885 


IC2944-122 


8.7 


0.54 


4.103 


HD184915 


13.8 


0.53 


3.791 


NGC869-1141 


9.7 


0.51 


3.824 



a The name of a cluster star has two parts, the cluster name 
and the WEBDA ID. 
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Fig. 1. — A plot of l^sini versus its fitting error in the low Vsini regime for our 2008 field 
B star sample. The solid lines indicate the mean value of A Vsini of the stars in each 20 
km s -1 bin. The dashed lines mark the positions that are larger than the AVsini values of 
90% of the stars in each bin. 
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Fig, 



A compa rison of the evolutionary tracks for non-rotat ing stellar models by 



Schaller et al.l ( 119921 ) (thick solid lines) and for rotating models by lEkstrom et al.l (120081 ) 
(thin lines: Q/Q ciit = 0.1; dashed lines: Q/Q crit = 0.5; dash-dotted lines: Q/Q crit = 0.9). 
The masses of the rotating models are 20, 9, and 3 M & from left to right, respectively. The 
masses (in M@) of non-rotating models are marked by the numbers at the bottom. 
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Fig. 3. — Cumulative distribution functions of projected rotational velocity for field (solid 
line) and cluster B stars (dotted line). 
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Fig. 4. — The logg po i ar distributions of single field (solid line) and cluster (dashed line) B 
stars. 
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Fig. 5. — Scatter plots of V sin i as a function of log g po iar for the field (left column) and cluster 
samples (right column) and for three mass ranges: 2 < M/M & < 4 (top), 4 < M/M & < 8 
(middle), and M/M Q > 8 (bottom). The solid horizontal lines indicate the mean Vsini 
of each bin containing six or more stars while dotted lines indicate the same for the bins 
containing fewer. Shaded areas illustrate the standard deviation of the mean for each bin. 
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Fig. 6. — The Vsini/V^i t histogram of all young B stars in our sample with logg po i ar > 4.15. 
Its polynomial fit is plotted as a thin solid line. The V cq /V cr i t distribution curve deconvolved 
from the polynomial fit is plotted as a thick solid line. 




Fig. 7. — Rotational velocity distributions (in the same format as Fig. [6]) for three subgroups 
of very young B stars with log 5r po iar > 4.15: 2 < M/M & < 4 (top panel), 4 < M/M Q < 8 
(middle panel), and M/M Q > 8 (bottom panel). 
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Fig. 8. — Scatter plots of (Vsini/V^ r i t , log g po iar) and rotational distributions (in the same 
format as Fig. [7]) for non-emission line stars. From top to bottom, each pair corresponds 
to mass ranges of 2.2 < M/M < 4.0, 4.0 < M/M < 5.4, 5.4 < M/M < 8.6, and 
M/M Q > 8.6. The shaded areas indicate the top 4% stars in the V cq /V cr i t distributions. 




Fig. 9. — Evolutionary changes in the EMMB models for M = 3M (top panel) and 9M Q 
(bottom panel). The solid curve lines are for models with the initial ratios of f2/f2 crit = 0.1, 
0.3, 0.5, 0.7, 0.8, 0.9, and 0.99 from bottom to top, respectively. The two vertical dotted lines 
in each panel mark the range in evolutionary state for the calculation of the mean velocity. 
The horizontal dotted line in each panel marks the rotational upper limit for normal B stars. 
The dashed lines in both panels are based on our case B models. 
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Fig. 10. — The V smi/V cvit vs. logg po i ar plots for B stars in mass range 2 < M/M & < 4 
(top panel) and in the range 7 < M/M & < 13 (bottom panel), using the same format as in 
Fig. El The thick curves a re theoretical p r edicti ons based on 3M (top panel) and 9 M & 
(bottom panel) models from lEkstrom et al.l (120081 ). The predictions of the case B model are 
plotted as thick dotted lines. The embedded plot in the bottom panel shows the rotation 
distribution curves for the YMS stars with 7 < M/M Q < 13 (same format as Fig. [6]). 



